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Abstract 



We propose a canonical tranformation approach to the effective interac- 
tion W e ff between two holes, based on the three-band Hubbard model but 
ready to include extra interactions as well. An effective two-body Hamiltonian 
can in principle be obtained including any kind of virtual intermediate states. 
We derive the closed-form analytic expression of the effective interaction in- 
cluding 4-body virtual states, describing the exchange of an electron-hole pair 
to all orders. The resulting integral equation, valid for the full plane, leads to 
a Cooper-like instability of the Fermi liquid. The two-hole bound states gen- 
eralize those reported earlier in cluster calculations by exact diagonalisation 
methods. 
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Typeset using REVTgX 



The three-band Hubbard model Hamiltonian has often been used to characterize the 
electronic properties of high-T,^ superconductors |1J as observed in electron spectroscopy 
experiments. The model is: 

H = H + W (1) 
where the independent hole hamiltonian reads, in the site representation 

Cu O n.n. 

where n.n. stands for nearest neighbors. The on-site repulsion Hamiltonian will be denoted 
by 

W = Ypim+m-, (3) 

i 

where Ui = Ud for a Cu site, Ui = U p for an Oxygen. We also considered first-neighbor 
0-0 hopping and off-site interactions 0, but they were dropped when it was clear that 
they were not really important. The point symmetry Group of the Cu-0 plane is C± v , and 
its characters are listed in Table I. We wish to show that holes are paired in a Cooper-like 
fashion in the ground state of this popular model. Preliminarly, however, we must introduce 
the W=0 pairs. Consider a determinantal 2-hole eigenstate of Hq built from degenerate 
Bloch states. Can we superpose such wave functions to obtain simultaneous eigenstates of 
H ,Ho and Wl For a general system, of course, we can't, but the symmetry of the problem 
enables us to achieve the result, and in addition to choose the eigenvalue W=0 of the on-site 
interaction. Indeed, omitting the band indices,let 

d[k] = \\k+, —k-\\ = c\ + cL k _\vac > (4) 

be a two-hole determinantal state derived from the Bloch eigenfunctions of H . We introduce 
the determinants Rd[k] = d[Rk] = d[kR], R G C^ v , and the projected states 

% W-^E X M (R) M[k] (5) 

where x {R) is the character of the operation R in the Irreducible Representation (Irrep) 
r). In the non-degenerate Irreps, the operations that produce opposite kn have the same 
character, and the corresponding projections lead to singlets. The explicit W=0 singlet pair 
states are: 

Xo 

$i B2 [A;, ri, r 2 ] = "^/|{ cos [M r i - r 2)} <P (K n) <p (k, r 2 ) 

- cos [kc^in - r 2 )] (k C4 , n) (k C4 ,r 2 ) 
- cos [k a (rx - r 2 )] <p (K, r x ) <p (k a , r 2 ) 
+ cos[k a ,(r 1 -r 2 ))(f>{k cr ',r 1 )(l)(k CT/ ,r 2 )} (6) 

and 
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$i A2 [k,r u r 2 ] 



Xo 
V2 



{cos [k(rt - r 2 )] r x ) (A;, r 2 ) 



+ cos [k CA {ri - r 2 )] (fc C4 , »"i) (&c 4J r 2) 
- cos [k a (ri - r 2 )\ (K, ri) (£v, r 2 ) 
- cos [Av(ri - r 2 )] (Av, r x ) (k a >,r 2 )}, 



(7) 



where xo is a singlet spin function and is the periodic part of the Bloch function. Using 
the explicit Bloch states we can verify that both vanish for r x = r 2 . Therefore they have 
the W=0 property. Turning now to the many-body problem, suppose the Cu-0 plane is in 
its ground state with chemical potential \i = E F and a couple of extra holes are added. We 
can consider the W=0 pairs as bare quasiparticles for which no repulsion barrier needs to 
be overcome for pairing. Do the dressed quasiparticles look like Cooper pairs? 

In Cooper theory ||, an effective interaction involving phonons is introduced via an 
approximate canonical transformation. Our problem is similar, except that the holes can 
exchange electron-hole pairs instead of phonons. Now the vacuum is the filled Fermi sphere 
and W connects the determinants d[k] to the 4-body (3 hole-1 electron) determinants that 
carry no quasi-momentum; they are of the form 



\a > 



(k' + q + k 2 ) + , k 2 _, -q_, -k'_ 



> 



(8) 



where k 2 is the electron state and pedices refer to the spin direction; those with opposite 
spin indices contribute similarly and yield a factor of 2 at the end. These are eigenstates of 
H : 



H \a) = E a | a) 
The interaction matrix element reads: 



(9) 



< 



(k' + q + k 2 ) + , k 2 -, ~k'_ \W\d\s\ >= 
5{q - s)U (q+ k' + k 2 , -k', s, k 2 ) 
-5 (k' — s)U (q + k' + k 2 , -q, s, k 2 ) 



(10) 



where the U matrix is computed on the Bloch basis. If we want to keep a close analogy 
with the Cooper theory, we temporarily consider W as a small perturbation and look for an 
approximate canonical transformation such that the new Hamiltonian 



H 



(11) 



contains no first-order terms. Here, H operates on the space of pairs, since the 3 holes- 1 
electron intermediate states have been decoupled by the canonical transformation. This is 
accomplished if 



W+[H ,A] =0 



(12) 



because then, to second order, 



H = H + -[W,A}. 



(13) 
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Assuming that the denominators do not vanish (more about that later) we obtain 

/ i a i \ (a\W\s) , n A . 

HA[s)= xl 1/ (14) 
Ei s — t, a 

We may write 

H — H + F + W e ff (15) 

where F is diagonal in the pair space, like Hq, and corresponds to the unlinked self-energy 
diagrams, while the effective interaction operator is W e ff. Like in Cooper theory, F will be 
dropped. We obtain 



2{p\F + W eff \s) 

1 1 

+ 



E p — E a E s — E a 



(16) 



Using the interaction matrix element (jlOf) , the product in (|16|) yields 4 terms; two are 



proportional to 5 (p — s) and belong to F, while the cross terms yield identical contributions 
to the effective interaction. 

After long algebra, we write the interaction between symmetry projected states (with s 
and p empty): 



where, explicitly, 



($ a [p}\W eff |$ a Is}) = J2 X {a) (R) (d [ P ]\W eff \Rd [s]) (17) 

R 



(d\p]\W eff \d[s])=2j2®(e(s + P + k)-E F ) 

k 

xll (s + p + k, —p, s, k) U (p, k, s + p + k, —s) 

1 



x 



e (s + p + k) - e (k) - e (s) + e (p) 

! , 

e(s +p + k) - e{k) -e{p) + e{s) 1 



(18) 



In the continuum limit, the integral must be understood as a principal part, because the 
denominators vanish in a domain of zero measure within the integration domain. However, 
this complication is an artifact of the perturbation approach, as will soon become clear. 
The approximate canonical transformation is open to question because it assumes that W 
be weak. Therefore, it is important to put the theory on a clearer and firmer basis, by 
removing the coupling to the a states to all orders. 

The transformation (0,0) corresponds to a second-order diagram for the two-hole am- 
plitude, and if one of the hole lines is closed on itself one gets the corresponding self-energy. 
The expansion of A can be continued systematically to produce the perturbation series. In- 
cluding all the diagrams of a generalized RPA would lead to something like the well-known 
FLEX approximation [|J whose implications for the superconductivity in the three-band 
Hubbard model have been explored recently in a series of papers || . A related self-consistent 
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and conserving T-matrix approximation has been proposed by Dahm and Tewordt |J for 
the excitation spectra in the 2D Hubbard model; we mention incidentally that recently 
diagrammatic methods have been successfully applied to the photoelectron spectra of the 
Cuprates in other contexts too, like the spin-fermion model 0. 

Realizing the key role of symmetry and W=0 pairs in this problem enhances intuition be- 
sides simplifying the perturbation formalism considerably Here, we wish to take advantage 
of the W=0 pairs to propose a non-perturbative approach to pairing based on a different way 
of performing the canonical transformation. We write the ground state wave function with 
two added holes as a superposition of two-body states (roman indices) and 4-body ones: 

l^o >= J2 a m\ m > +XX|a > • (19) 

m a 

Schrodinger's equation then yields 

(E m — E )a m + ^a m iV mtm i + y~"b a W m , a = (20) 

m! a 



(E a - Eo) b a + 5> m ,W a , m , = (21) 

rn' 

where Eq is the ground state energy. V m i, m vanishes for W=0 pairs in our model; however 
we keep it for generality, since it allows to introduce the effect of phonons |||9| and any other 
indirect interaction [|IO|-|i"2|j that we are not considering. Solving for b a and substituting in 
(p0|) we exactly decouple the 4-body states. The eigenvalue problem is now 

(Eq — E = J2 a m'{Vm, m ' + (m\S[E }\m')}, (22) 

m' 

where 

, ^< m\W\a >< a\W\m' > , , 

(m\S [E ] \m') = £ !— L- -J— ! . (23) 

a E - E a 

This is of the form of a Schrodinger equation with eigenvalue Eq for pairs with an effective 
interaction V + S. Then we interpret a m as the wave function of the dressed pair, which 
is acted upon by an effective hamiltonian H. This is the canonical transformation we were 
looking for. However, the scattering operator S is of the form S = W e ff + F, where W e ff 
is the effective interaction between dressed holes, while F is a forward scattering operator, 
diagonal in the pair indices m ,m' which accounts for the self-energy corrections of the one- 
body propagators: it is evident that it just redefines the dispersion law E m , and, essentially, 
renormalizes the chemical potential. Therefore F must be dropped, as in Cooper theory 
and above. Therefore the effective Schrodinger equation for the pair reads 

(H + V + W eff )\a>= E \a> (24) 

and we are interested in the possibility that — ^Fp — |A|, with a positive binding energy 
|A| of the pair. The V interaction just adds to W e ff, and this feature allows to include 
in our model the effects of other pairing mechanisms, like off-site interactions, inter-planar 
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coupling and phonons. Again, the product in the numerator of fl23|) yields 4 terms; two are 
proportional to 5(p — s) and belong to F, while the cross terms yield identical contributions 
to W e ff. Hence we obtain the effective interaction between W=0 pairs: 

occ 

<$„ [p}\ W eff \% [s]) = 4 £ X {v) (R) E ( £ (Rs+p + k)- E F ) 

R&Ca v k 

U (Rs + p + k, —p, Rs, k) U (p, k, Rs + p + k, —Rs) 
X £ (Rs + p + k) - e (k) + e (s) + e (p) - E ( ' 

The sum is over occupied k with empty Rs + p + k. There are no vanishing denominators 
in the new expression, if E is off the continuum. The p and s indices run over 1/8 of the 
Brillouin Zone. We denote such a set of empty states e/8, and cast the result in the form of 
a (Cooper-like) Schrodinger equation 

e/8 

2e (k) a(k) + ^2 w eff (K k') a (k') = E a (k) (26) 
k' 

for a self-consistent calculation of E (since W e ff depends on the solution). Let N c be the 
number of cells in the crystal. The U matrix elements scale as N^ 1 and therefore W e ff 
scales in the same way. For an infinite system, Nc oo , this is a well defined integral 
equation. In principle, we can do better. By a canonical transformation one can obtain 
an effective Hamiltonian which describes the propagation of a pair of dressed holes, and 
includes all many-body effects. The exact many-body ground state with two added holes 
may be expanded in terms of excitations over the vacuum (the non-interacting Fermi sphere) 
by a configuration interaction: 

|^o >= E fl ™l m > > + E^l/ 3 > +•••• ( 27 ) 

here m runs over pair states, a over 4-body states (2 holes and 1 e-h pair), {3 over 6-body 
ones (2 holes and 2 e-h pairs), and so on. To set up the Schrodinger equation, we consider 
the effects of the operators on the terms of >• We write: 

H \m >= E m \m >, H \a >= E a \a >, ... (28) 

and since W can create or destroy up to 2 e-h pairs, 

W\m >= Y^ V m',m\m' > +J2\ a > W *,m 
m' a 

+ El/5 > Wp, m , (29) 

P 

and the like. The Schrodinger equation now gives an infinite system for the amplitudes of 2n- 



body states; however we can show []13| that these amplitudes can be successively decoupled, 
producing a renormalization of W matrix elements and energy eigenvalues E a , Ep and so on. 
In principle, the method applies to all the higher order interactions, and we can recast our 
problem as if only 2 and 4-body states existed. If one calculates W e ff neglecting 6-body and 
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higher excitations, at least the structure of the solution is exact when expressed in terms of 
renormalized matrix elements. 

The Hubbard model, with V=0, leads to pairing, and the mechanism was first discovered 
in cluster studies @. Now we can better understand those results, since the above theory- 
applies not only to the plane but also to clusters, provided that they are fully symmetric and 
allow W=0 pairs. The symmetry of the cluster is essential, because only fully symmetric 
clusters allow such solutions. The planar lattice structure is also essential, because no W=0 
pairs occur in 3D or in a continuous model. We studied PJT^-[T6| the fully symmetric clusters 



with up to 21 atoms by exact diagonalisation. The main difference is that in the clusters the 
symmetries of W=0 pairs were found [|17|,|T8f to be 1 B 2 and X A\. The reason for having 1 A 1 



instead of 1 A 2 is a twofold size effect. On one hand/Ai pairs have the W=0 property only 
in the small clusters, having the topology of a cross, and belonging to the S4 Group, but do 
not generalize as such to the full plane, when the symmetry is lowered to C± v ; on the other 
hand, those small clusters admit no solutions of X A 2 symmetry. To monitor pairing, initially 
we used || a definition of A taken from earlier cluster calculations [|T9| , |2"0|| (where, however, 



only unsymmetric clusters were considered, and this mechanism could not operate): 

A = E(N + 2) + E(N)-2E(N + 1), (30) 

where E(N) is the ground state energy of the cluster with N holes, as obtained by exact 
diagonalisation. Pairing, that is, A < 0, was found when (and only when) the least bound 
holes formed a W=0 pair. In small clusters this conditions the occupation number, because 
the holes must partly occupy degenerate (x,y) orbitals, while in the full plane the W=0 
pairs ([|,|7]) exist at the Fermi level for any filling. The second-order approximation A^ 2 ^ was 
obtained by ground state energy diagrams (modified for degenerate ground states, when 
appropriate); the resulting expression |L7| is clearly just the second-order approximation to 



the diagonal terms of Q2"5|). Then, we computed to second-order the two- hole amplitude for 
holes of opposite spins in the degenerate (x,y) orbitals. We demonstrated that this produced 



A( 2 ) 
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an effective interaction, which pushes down the singlet and up the triplet by 
agreement between the second-order calculation and the numerical exact diagonalisation 
results supported the interpretation. Thus, we have shown that in the symmetric Cu-0 
clusters genuine pairing takes place, due to an effective interaction which is attractive for 
singlets and repulsive for triplets. The cluster calculations [I3-IB| showed that W=0 pairs 



are the "bare" quasiparticles that, when "dressed", become a bound state. That approach 
is inherently limited by the small size of solvable clusters, but allows a very explicit display 
of paired hole properties, that even show superconducting flux-quantization |fT?| , Pj}] . 

The equations (^,^) allow to extend the study of pairs to the full plane. The integrand 
is very discontinuous because of Umklapp terms and of the limitations to occupied or empty 
states; moreover, its angular dependences are involved. The numerical solution is far from 
trivial, and the method that we developed will be presented elsewhere, with the numerical 
results. The main finding is the instability of the Fermi liquid in the model at hand. For 
both l B 2 and X A 2 we find pairing with a doping-dependent A. The three-band Hubbard 
model might be too idealized to allow a detailed comparison with experiments; however we 
stress that contributions from phonons and other mechanisms can be included by a non zero 
V. The approach presented here is more general than the model we are using, and can be 
applied to realistic Hamiltonians. 
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TABLES 



C4U 


E 


c 2 


2C 4 


2a 


2a' 




Ai 


1 


1 


1 


1 


1 




A 2 


1 


1 


1 


-1 


-1 


R z 


Bi 


1 


1 


-1 


1 


-1 


x 2 — y 2 


B 2 


1 


1 


-1 


-1 


1 


xy 


E 


2 


-2 














TABLE I. The Character Table of the C Av Group 
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